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Motivated by the recent observational and theoretical evidence that long Gamma-Ray Bursts 
(GRBs) are likely associated with low metallicity, rapidly rotating massive stars, we exam- 
ine the cosmological star formation rate (SFR) below a critical metallicity Zcrit ~ Zq/10 — 
I Zq/5, to estimate the event rate of high-redshift long ORB progenitors. To this purpose, we 

exploit a galaxy formation scenario already successfully tested on a wealth of observational 
data on (proto)spheroids, Lyman break galaxies, Lya emitters, submm galaxies, quasars, and 
local early-type galaxies. We find that the predicted rate of long GRBs amounts to about 
300 events yr^^ sr^^, of which about 30 per cent occur at z > 6. Correspondingly, the GRB 
QQ I number counts well agree with the bright SWIFT data, without the need for an intrinsic lu- 

. minosity evolution. Moreover, the above framework enables us to predict properties of the 

■ GRB host galaxies. Most GRBs are associated with low mass galaxy halos A/h ^ 10^^ Mq, 
, and effectively trace the formation of small galaxies in such halos. The hosts are young, 

■ with age smaller than 5 x 10^ yr, gas rich, but poorly extincted (Ay ^ 0.1) because of their 
OO ! chemical immaturity; this also implies high specific SFR and quite extreme a-enhancement. 

■ Only the minority of hosts residing in large halos with Afn ^ 10^^ Mq have larger extinction 

(Ay ~ 0.7 — 1), SFRs exceeding 100 A/© yr~^ and can be detected at submm wavelengths. 
Most of the hosts have UV magnitudes in the range — 20 ^ Afiaso ^ — 16, and Lya luminosity 

in the range 2 x 10*° < Lhya ^ 2 x 10*^ erg s~^. GRB hosts are thus tracing the faint end of 
' the luminosity function of Lyman break galaxies and Lya emitters. Finally, our results imply 

that the population of 'dark' GRBs occur mostly in faint hosts at high redshift, rather than in 
dusty hosts at low redshift. 
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1 INTRODUCTION 

The spectroscopic detection of the energetic supernova SN 2003dh 
coincident with GRB 030329 (Hjorth et al. 2003; Stanek et al. 
2003) has firmly established that - at least some - long GRBs ac- 
company the core collapse of massive stars, as it was first suggested 
by the spatial and temporal coincidence of GRB 980425 and SN 
1998bw (Galamaet al. 1998). Such spectroscopic signatures of su- 
pemovae (SNae) associated with GRBs have been detected in a 
handful of cases during the last years, e.g. GRB 031203/SN 20031w 
(Malesani et al. 2004), GRB 02121 1/ SN 20021t (Delia Valle et al. 
2003), GRB 050525A/SN 2005nc (Delia Valle et al. 2006), and 
the recent case of GRB 060218/SN 2006aj (Campana et al. 2006, 
Modjaz et al. 2006). The possibility that core-collapse SNae are 
progenitors of GRBs has been further supported by the detection of 



re-brightening in the late-time afterglow light curves, interpreted 
as a contribution of accompanying SNae (Bloom et al. 1999, Zeh 
et al. 2004, Castro-Tirado & Gorosabel 1999) and the localizations 
of afterglows in star forming regions (e.g. Djorgovski et al. 2001, 
Bloom et al. 2002, Fynbo et al. 2000, Fruchter et al. 1999, Holland 
& Hjorth 1999). 

Indeed, the currently mostly favored scenario for the origin of 
long GRBs involves the collapse of a massive Wolf-Rayet star en- 
dowed with rotation (Woosley 1993, MacFadyen & Woosley 1999, 
Woosley & Heger 2004). Recently, Yoon & Langer (2005) con- 
sidered the evolution of massive, magnetized stars where rapid ro- 
tation induces an almost chemically homogeneous evolution, and 
found that the requirements of this collapsar model are satisfied if 
the metallicity is sufficiently small, namely less than 0.1 Zq (see 
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Woosley & Heger 2006; Yoon et al. 2006). This is broadly consis- 
tent with the estimates of metallicities of long GRB hosts, which 
yielded preferentially subsolar - down to 10~^ Zq - values (e.g. 
Gorosabel et al. 2005; Chen et al. 2005; Starling et al. 2005). 

Clearly this scenario has key implications not only on the 
physics of the event and on the evolution of massive stars at low 
metallicities, but also on the event rate and redshift distribution of 
long GRBs compared with that of SNae. The cosmological conse- 
quences of a metallicty threshold can be explored by considering an 
average cosmic metallicity evolution, as in the work by Langer & 
Norman (2006). In the present paper we instead explore the effects 
of a metallicity threshold following the star formation and chemical 
evolution of individual galaxies. 

In order to derive the rates of progenitors and the character- 
istics of their host galaxies at high redshift, the star formation his- 
tory and the chemical evolution for a large range of galaxy mass 
and virialization redshift must be computed. Baryon condensation 
in cold gas, stars and in a central massive black hole (BH) within 
galaxy Dark Matter (DM) halos is a quite complex outcome of a 
number of physical processes (including shock waves, radiative and 
shock heating, viscosity, radiative cooling; star formation, BH ac- 
cretion, gas inflow and outflow) largely affecting each other (see 
Granato et al. 2001, 2004; Croton et al. 2006; De Lucia et al. 2006). 
When treating these processes within a self-consistent cosmologi- 
cal framework of galaxy formation, most of the complexity is re- 
lated to the different scales involved. 

For the sake of definiteness, here we adopt as a reference 
the galaxy formation scenario developed by Granato et al. (2004), 
which consistently accounts for the coevolution of spheroidal 
galaxies and their nuclear activity and intrinsically follows in time 
the gas content, star formation rate (SFR) and metallicity evolution 
for each galaxy mass. This enables us to investigate in detail the 
effect of a metallicity threshold on the properties not only of the 
GRB population, but also of their host galaxies. 

The outline of this paper is the following. In § 2 we briefly 
review the conceptual issues of the adopted galaxy formation sce- 
nario and explain how the GRB progenitor rate has been estimated. 
Our results are presented in § 3. In § 4 we discuss our findings, by 
comparing them with observational results and with previous stud- 
ies. In § 5 we summarize our conclusions. Throughout the paper a 
flat cosmological model with matter density parameter Om = 0.27 
and Hubble constant Ho = 72 km s~^ Mpc~^ is adopted. 



2 MODELLING 

2.1 Overview of the galaxy formation scenario 

Long GRBs at high redshift (z > 1) have progenitors which formed 
at least 8 Gyr ago. Their coeval stellar populations are as old as 
the populations of spheroidal galaxies and spiral bulges and older 
than the populations in present galaxy discs (see Renzini 2006): 
thus high redshift long GRBs trace the formation history of the 
oldest stellar populations. Low redshift GRBs conversely trace the 
star formation in small mass irregular/interacting galaxies with low 
SFR (Fruchter et al. 2006; Wainwright et al. 2007). Since we focus 
on high-z GRBs, in the following we will neglect the contribution 
of star formation in discs of present-day spiral galaxies. 

We exploit the physical model elaborated by Granato et al. 
(2001, 2004), which follows the evolution of baryons within pro- 
togalactic spheroids taking into account the effects of the energy 
fed back to the intra-galactic gas by SN explosions and by accre- 
tion onto the nuclear supermassive BH (see also De Lucia et al. 



2006; Croton et al. 2006). The model envisages that during or soon 
after the formation of the host DM halo, the baryons falling into 
the newly created potential well are shock-heated to the virial tem- 
perature. The hot gas is (moderately) clumpy and cools fast espe- 
cially in the denser central regions, yielding a strong burst of star 
formation. Star formation also promotes the storage of the cooled 
gas into a reservoir around the central seed BH, eventually leading 
to accretion onto it (see Kawakatu & Umemura 2002). The ensu- 
ing SN explosions and the nuclear activity feed energy back to the 
baryons, and regulate the ongoing SFR and BH growth. These mu- 
tual energy feedbacks actually reverse the formation sequence of 
the baryonic component of galaxies compared to that of DM halos: 
the star formation and the buildup of central BHs are completed 
more rapidly in the more massive haloes, thus accounting for the 
phenomenon now commonly referred to as downsizing (e.g. Cowie 
et al. 1996; Kodama et al. 2004; Glazebrook et al. 2004). 

In Appendix A we present a simplified version of basic equa- 
tions of the model, which allow to derive analytical solutions for 
the time evolution of SFR, mass in stars and metallicity, the quanti- 
ties relevant to this work. These analytical functions are very good 
approximations of the more complex system of equations numeri- 
cally solved in Granato et al. (2004; see for details Lapi et al. 2006 
and Mao et al. 2007). Because of their fundamental character, these 
equations catch the basic aspects of the physical processes ruling 
star formation in protogalaxies at high redshift. The adopted Initial 
Mass Function (IMF) is a double power-law with slope 1.25 from 
120 Mq to 1 Mq and 0.4 from 1 Mq down to 0.1 M© (Romano 
et al. 2002), which is quite similar to that proposed by Chabrier 
(2005). 

In Fig. 1 we present the results for DM halos virialized at 
z = 6 and endowed with mass ranging from 10"' to 10^^ Mq. 
For masses Mh > 10^^ Mq the SFR increases almost linearly with 
galaxy age in the initial stages, and then it is suddenly halted by the 
energy feedback from the quasar after a few 10* yr. On the contrary, 
for Mh 10^^ Mq it is first almost constant and then slowly de- 
clines due to gas exhaustion. In fact, for small masses the effect of 
SNae feedback regulates the star formation, while the BH is rather 
small there and the SFR can proceed for a much longer time (see 
Eq. [16] in Shankar et al. 2006). 

The mass cycled through stars at any time is easily obtained 
by integrating the SFR (see Fig. 1). Note that the quantities plotted 
in Fig. 1 refer to redshift z = 6; however, both SFR and mass in 
stars scale approximately as (1 -I- z)^^^ for a given halo mass. 

From the figure it is apparent that the chemical enrichment 
of the cold gas component is very rapid; e.g. with the adopted 
IMF the gas attains 1/100 and 1/10 the solar abundance in about 
1.2 X lO'^ and 5 x lO'^ yr respectively, almost independently of 
the halo mass and redshift. A Salpeter IMF yields timescales about 
twice as long. This rapid enrichment is due to the fast evolution 
and formation timescales (t*, see Appendix A) for the massive 
stars (greater than 10 Mq) relative to the timescale for the infall of 
the diffuse mediimi with primordial composition, which dilutes the 
metallicity of the cold star-forming gas. This behavior shows that 
a possible change of the IMF at metallicity lower than a threshold 
around 3 — 5 x 10^'' Zq (see Bromm et al. 2001; Schneider et al. 
2006) is not critical for the issue related to long GRBs. In other 
words, if long GRBs are associated with low metallicity environ- 
ments, the quick enrichment implies that the most relevant epoch 
for GRB progenitors comes soon after halo virialization. 

This galaxy formation model neglects spatial resolution and 
assimies instantaneous mixing, i.e. it averages both SFR and chem- 
ical composition over the entire mass of cold gas. We stress that the 
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Figure 1. SFR, stellar mass, specific SFR, average cold gas metallicity, cold 
gas mass, and extinction at 1350 A (from top to bottom) as a function of 
the galactic age for halos of masses 10^" Mq {solid lines), 10^^ Mq (dot- 
dashed hnes), 10^^ Mq {dashed lines), 10^^ Mq (c/offerf lines), virialized 
at redshift z = 6. The vertical lines mark the epoch when the cold gas 
metallicity attains the critical thresholds Zq/10 {green) and Zq /5 (red). 



cold gas is only a small fraction of the overall baryons associated 
with the galaxy halo (cf. Appendix A), and thus the averaging con- 
cerns the mass/volume of the protogalaxy wherein star formation is 
occurring and not the overall DM halo mass/volume. On the other 
hand, metallicity gradients have been observed in local spheroidal 
galaxies, A[Z/H]/Alogr ^ -0.25 (e.g. Annibah et al. 2007, 
Sanchez-Blanquez et al. 2007). This implies variations of about a 
factor two in metal abundance within a radius enclosing most of 
the galaxy mass. As we shall see such a factor is not crucial to the 
conclusions of this paper. 

The model here adopted has the very valuable asset that 
successfully fits a wealth of observational data and constraints 



regarding protospheroids, Lyman break galaxies and Lya emit- 
ters, submm-selected galaxies, quasars, EROS and local early-type 
galaxies; for a detailed comparison with the observational data, we 
defer the interested reader to the papers by Granato et al. (2001, 
2004, 2006), Cirasuolo et al. (2005), Silva et al. (2005), Lapi et al. 
(2006, in particular their Table 2), Mao et al. (2007); furthermore, 
clustering properties of submm galaxies have been extensively dis- 
cussed in the context of our model by Negrello et al. (2007). 

2.2 Estimate of long GRB progenitor rates 

The key ingredients provided by the galaxy formation model re- 
lated to this work are the SFR A/* (t) and the average gas metallic- 
ity Z{t) as a function of the age t for an individual galaxy within 
a given halo mass Mh virialized at redshift z (see Appendix A for 
handy approximations). Notice that Z refers to the average metal- 
licity of the gaseous component from which new stars form. 

The cosmic SFR per unit volume at redshift z contributed by 
objects with average gas metallicity Z < Zcvit is thus given by 



Z<Z„it 



dM*dMH 



d^iVsT d 
dMndi^ dlnM^ 



(1) 

where d^ Ai'sT /dA/n dfz are the formation rates of DM halos at 
cosmic time computed using the Sheth & Tormen (1999) mass 
function, and 7^^^* is the time the galaxy spends at SFR higher 

than M* and metallicity lower than Zcrit • 

It is plausible that the core-collapse, believed to give rise to 
the formation of a BH and a GRB event, is related to the final 
phases of the evolution of stars with masses greater than 12 Mq 
and with rapidly rotating cores (Yoon & Langer 2005; Woosley & 
Heger 2006; Yoon et al. 2006). From a theoretical point of view a 
major issue in this respect is the large mass loss, that would entail 
also large angular momentum loss. This problem is possibly alle- 
viated in stars of metallicity below a critical threshold lower than 
1/5 — 1/3 Zq and high initial spin rate (Woosley & Heger 2006; 
Yoon et al. 2006). Though observations are still scanty in support- 
ing this expectation, we will consider low metallicity and high ro- 
tation velocity as necessary and sufficient conditions for a massive 
star to be considered as a GRB progenitor. 

The GRB progenitor rate will thus be based on the above pre- 
liminary theoretical estimates which suggest that the expected frac- 
tion of GRB progenitors with respect to the number of massive stars 
(8 Mq < m* < 100 Mq) is /prog > 2 - 3 per cent for metallicity 
Zcrit < ^0/5 (see Woosley & Heger 2006; Yoon et al. 2006; see 
also Bissaldi et al. 2007). 

The absolute GRB progenitors rate per unit volume can be 
then expressed as 

/prog WSN STTZ[z)z<z, 

nlOO 



7^o 



(2) 

dm* is the num- 



k-prog(2) 

where nsN = /g""™ (/>(m*)dm*/ J^"^' mi,(p[mi, 
ber of massive stars ending in SNae per unit mass of formed stars; 
for the adopted IMF nsN ~ 1-4 x 10"^ (it halves for 

a Salpeter IMF). Thus assuming /prog ~ 0.02, the number of 
GRB progenitors per unit mass of formed stars amounts to about 
3.6 X 10"* M-'^ for Zcrit < ^o/5. 

2.3 Long GRB number counts 

The number of GRBs at redshift greater than z is 

/ Tiprosiz') AV 



7^: 



GRB 



(> Z) ^ /be 



dz' 



1 + z' dz' 



(3) 
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where V is the cosmological volume per unit solid angle and the 
factor (1 + z)~^ accounts for time dilation effects due to red- 
shift. Moreover, recall that GRBs are believed to be anisotropic 
phenomena owed to the flow collimation and/or relativistic beam- 
ing effects. Thus only a fraction /beam of the estimated progen- 
itors prompt emission would point toward us within an opening 
angle 2 9 and could then be observed as a GRB. For the events 
for which a (jet) opening angle have been estimatetfB assuming 
that an afterglow light curve break owed to a jetted structure with 
a half angle 9 (Ghirlanda et al. 2007), the median value ^ ?a 6 
deg has been inferred (see also Guetta et al 2005); correspond- 



ingly^ /be 



5.5 X 10 (^/6deg) . The expected number 



of long GRBs per unit mass of formed stars within metal poor 
environments (Z ^Zcvit) is thus k ~ 1.5 x 10"'' (nsN/1.4 x 
10-') (/prog/0.02) (9/6 deg)2 

The long GRB number counts at flux limit S\\m and redshift 
> 2 is given as 



^grb(> z)s>Sii 



/b 



dz — — — I P{L)dL ; 



1 



(4) 

in this expression P{L) is the equivalent isotropic luminosity dis- 
tribution (see § 3.2 for details), and is the luminosity corre- 
sponding to the limiting flux, given by 



_ 47r dljz) 



(5) 



in terms of the luminosity distance dL{z) and of the 7f -correction 
K(z); note that the latter quantity depends on the GRB spectra, see 
§ 3.2 for details. 



3 RESULTS 

3.1 Long GRB progenitor rates and redshift distribution 

If progenitors of long GRBs are metal poor massive stars, the metal 
abundance Z of the cold gas wherein stars form is a crucial phys- 
ical parameter. Fig. 2 shows the cosmological SFR for the over- 
all galaxy population and for systems with Z^Zcvit ~ 0.1 and 
0.2 Zq as function of redshift. The model prediction can be fitted 
with the approximate formula 



, / STTZ(z 
log ( 



Moyr-^Mpc" 



^a + b{z-Zmax) , 
3.8 (-3.2) X 10"' and Zn 



(6) 



where a = -1.4 (-3.2), 6 = - 
4.5 (6.5) for Zcm = (^©/lO). 

The predicted overall SFR reproduces fairly well the 
extinction-corrected data at 2 < z < 6, with reasonable correction 
for dust extinction. At z < 1 the observed SFR is underestimated 
due to the fact that intrinsically our model does not account for the 
contribution of star formation in discs of present spiral galaxies. On 
the other hand, only a minor fraction (less than 20 per cent) of long 
GRBs with redshift determination (though it is difficult to quantify 
the observational biases which certainly affect this percentage) are 
located at 2: < 1 , while we are interested in the bulk of the (high 
redshift) GRB population (cf. § 2.1). 

As expected the SFR in GRB hosts decreases with Zait, 
but the effect is differential with redshift. Since the timescale 

^ These estimates maybe clearly affected by selection effects, though not 
easy to quantify. 
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Figure 2. Cosmological SFR (left y-axis) and rate of long GRB progeni- 
tors (right y-axis) as a function of redshift, computed with no threshold on 
metallicity {blue line), and with thresholds at Zq/^ (red line) and Zq/10 
(green line). Data are from Schiminovich et al. (2005; circles), Steidel et al. 
(1999; squares), and Bouwens et al. (2006; stars). The shaded area illus- 
trates the uncertainties due to the extinction correction; the data have been 
rescaled down by a factor of about 2 to account for the adopted IMF (see 
text). In the hatched region the SFR is dominated by discs of present spiral 
galaxies, not treated in our model. 




Figure 3. Normalized redshift distribution of GRB progenitors, computed 
with no thi'eshold on metallicity (blue line), and with thresholds at Zq / 5 
(red line) and Zq/10 (green line). 



tcrit ^ 5 X 10*^ yr to reach the critical metallicity threshold is al- 
most independent of mass and redshift (cf. Fig. 1), its ratio to the 
cosmic time increases toward high z for all halos. The fraction of 
overall cosmic star formation occurring in metal poor protogalaxies 
raises with redshift, and so does the predicted rate of GRB progeni- 
tors. The net effect is that the expected redshift distribution of long 
GRB progenitors peaks to a redshift significantly higher than that of 
the cosmic SFR. As a consequence the higher the redshift, the more 
directly GRBs mirror the cosmic SFR. The fall off of star formation 
beyond z > 10 in the galaxy formation model is a consequence of 
the effect of the SN feedback; at increasing 2 the decrease of star 
formation oc is not balanced by the rise in the number of 

virialized halos. 
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Figure 4. GRB counts in the 15 — 150 keV band, computed with no thresh- 
old on metallicity (blue line), and with thresholds at Zq/^ (red line) and 
Zq/10 (green line). Data (filled dots) are from the 2-yr SW/Fr catalog. 
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z 

Figure 5. The redshift distribution of GRBs at the hmiting flux [15-150 
keV] of 0.01 (dotted line), 0.1 (dashed line), and 1 ph s^^ cm~^ (solid 
line), computed with a threshold on metallicity at Zq/W. 



The long GRB progenitor rate shown in Fig. 2 can be inte- 
grated over cosmic time to obtain the normalized redshift distribu- 
tion reproduced Fig. 3. For the case ^crit = Zq/1Q, the fraction of 
progenitors is around 60 per cent at z > 4 and around 30 per cent 
at z > 6. 

If the threshold Zc-cit were to refer to stellar rather than ISM 
metallicity, the time within which GRB progenitors can be pro- 
duced in any galaxy halo would be overestimated by a factor around 
2. This would in turn double the duty cycle and, hence, the num- 
ber of progenitors. The difference in the predictions worsens with 
higher Zcrit threshold, to the point where the threshold might never 
be reached in small galaxies, which then can host GRB progenitors 
over the whole Hubble time. Again the effect is differential with 
redshift, with a larger increase in the progenitor number at higher 
z. 

3.2 Long GRB counts with SWIFT 

In order to assess whether a scenario in which single metal poor, 
rotating massive stars are the GRB progenitors is tenable within the 
adopted galaxy model framework, we estimated the predicted long 
GRB counts which should have been detected by SWIFT (Gehrels 
et al. 2004). 

As discussed in § 2.2 and § 2.3 we consider a number of ob- 
served GRBs per unit mass k ~ 1.9 x 10~^ Af© ^ ™d 3 metallicity 
threshold Zcnt ~ 1/10 - 1/5 Zq. 

In order to estimate the observable number of GRBs as func- 
tion of limiting flux, we assume that the prompt GRB luminosity 
distribution can be characterized, independently of redshift, as 

P(L) cxL-*e-^<^/-^ , (7) 

with a low luminosity cutoff at Lc = 3 x 10^^ erg s~^, and slope 
5 = 2. This parametrization is consistent with those constrained by 
Daigne et al. (2006) and Guetta et al. (2005). As for the spectrum, 
we adopted - as commonly done - a typical Band representation, 
with low and high energy slopes q = — 1 and (3 = —2.25, re- 
spectively. The peak energy is considered to follow a log-normal 
distribution (Preece et al. 2000) with mean log i5pcak,o=2.74 and 
dispersion a = 0.3 dex (see Daigne et al. 2006). Although such 
spectral parameters have been estimated from bright BATSE GRBs 
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Figure 6. The redshift distribution of GRBs computed with a threshold on 
metallicity at Zq/10 (solid line). The other curves illustrate the contribu- 
tions from host galaxies with different halo masses: 10^ — 10^" Mq (triple- 
dot-dashed line), 10^° - 10^1 Mq (dot-dashed line), lO" - 10^^ j^^ 
(dashed line), lO" - lO^^ Mq (dotted line). 

(Preece et al. 2000), this provides the simplest hypothesis which - 
given the relatively small effect owed to the iC-correction - appears 
adequate for the consistency check we intended to perform. 

The predicted number counts are shown in Fig. 4 together with 
the actual GRB SVWFr counts. A meaningful comparison can only 
consider 'bright' SWIFT GRBs, namely above a photon flux of 1 
ph s~^ cm~'^, as below this level some degree of incompleteness is 
expected (Band 2006). 

As it is apparent from Fig. 4 a reasonable agreement with 
the SWIFT counts can be obtained after the above assumptions. It 
should be pointed out that we did not perform any fitting optimiza- 
tion, but simply compared the counts predicted by the model with 
the data, under the simplified and commonly adopted assumptions 
on the luminosity and spectral energy GRB distributions outlined 
above. 

A very interesting aspect concerns the fact that the counts can 
be reasonably reproduced without requiring any GRB prompt lu- 
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Figure 7. Extincted magnitude at 1350 A and Lya luminosity of the host 
galaxies of the GRB progenitors as a function of galactic age. Lines as in 
Fig. 1. The first and third panels refer to redshift 3, the second and forth 
ones to redshift 6. The shaded areas illustrate the ranges sampled in the 
observed luminosity functions. 



minosity evolution (e.g. Daigne et al. 2006). This is a natural con- 
sequence of the fact that our redshift progenitor distribution intrin- 
sically peaks at redshift higher than the cosmic star formation. In- 
deed, if we were to reproduce the number counts with no metallic- 
ity threshold, the number of GRBs per unit mass of formed stars 
(fc) would have to be a factor about 10 lower, which would imply 
that the bulk of the progenitors would be located at lower redshift. 

The predicted redshift distribution of GRB events for different 
flux limits is presented in Fig. 5. This shows that already at fluxes 
5* > 1 ph s~^ cm~^ we expect that 10 per cent of GRBs occur at 
2 > 6. Of course such a fraction increases with decreasing limit- 
ing flux, reaching around 30 per cent for S > 10^^ ph s^^ cm^^; 
since this flux limit is below the present detectable level, the upper 



curve of Fig. 5 represents the redshift distribution of all currently 
observable GRBs. 

More specifically, the model predicts that in two years 13 
GRBs at 2;>6 and flux limit 1 ph s^^ cm^^ should have been 
detected by SWIFT. So far only 1 GRB has robust redshift esti- 
mate at 2 > 6. On the other hand, while it is difficult to quantify 
the efficiency of the redshift determination, especially for very high 
redshift events, this could well be less than 10 per cent. 



3.3 GRB host galaxies 

The exploited galaxy formation model also enables us to directly 
predict the properties of the GRB host galaxies (SFR, magnitude, 
stellar mass, average metallicity, extinction) as function of the halo 
mass and age as well as the luminosity function (LF) and the cor- 
responding SFR distribution of the overall population. 

As discussed in § 2. 1, during the evolution of individual galax- 
ies, the metallicity of the star-forming gas attains the threshold 
Zcrit = 0.1 — 0.2 Zq quite rapidly, within tcrit ~ 5 x 10*^, and 
this timescale is basically independent of the host halo mass. 

Interestingly, since this is much shorter than the halo virializa- 
tion timescale, tait <C tn, the GRB rates for fixed halo mass di- 
rectly trace the cosmic rate of halo virialization. As a consequence, 
the GRB hosts mostly reside within galaxy halos in the mass range 
10^ < Mh < 10" Mq, and only at 2 < 4 the fraction of hosts in 
massive halos Mu ^ 10^^ Mq exceeds 10 per cent. The relative 
number of GRBs for different halo masses as function of redshift is 
shown in Fig. 6. 

The predicted evolution of the SFR and stellar mass as func- 
tions of galactic age and halo mass are reported in Fig. 1. For halo 
masses in the range 10^° - 10^^ Mq and virialized at redshift z^6 
the expected SFR spans 0.3 — 10^ Mq yr~^, and the corresponding 
stellar masses cover the interval 10^ < < 2 x 10^° Mr.-,, with 

both quantities scaling as (1 + 2)^''^ at a given halo mass. The 
shortness of t„it yields high specific SFRs, A'L/M^, > 2 x 10"**, 
almost independently of the virialization redshift. 

Now we turn to consider more directly observable properties 
of GRB hosts, such as the average ultraviolet extinction ^1350 and 
the corresponding extincted magnitude M1350 at 1350 A, again 
as function of galactic age. We recall that M1350 ~ —18.6 — 
2.5 log(M*/MQ yr"^) + ^1350 (see Appendix A). Fig. 1 reports 
^1350, which scales with redshift as (1 + 2)^''^, for various halo 
masses virialized at redshift 2 = 6. A key prediction of the model 
is apparent, namely that most GRB hosts are poorly extincted sys- 
tems. Since the less massive hosts (Mh ^ 10^^ Mq) largely out- 
number the most massive ones, the typical j4i35o ranges between 
0.01 — 0.3 mag, corresponding to Ay < 0.1. Larger dust extinc- 
tion is predicted only for more massive hosts Mh > 10" Mq, 
which exceed a few per cent of the total GRB hosts only at 2 < 6 
(cf. Fig. 1). 

The AB absolute magnitude at 1350 A, M1350, is reported in 
Fig. 7 for various halo masses and redshifts 2 = 3 and 2 = 6. 
The shaded area indicates the absolute magnitude range where the 
UV high-2 LF is currently well sampled. It is clear that the UV 
luminosity is practically not affected by dust, since the predicted 
extinction is low for galaxy ages less than tcrit ■ 

Fig. 7 also shows the age dependence of the expected Lya 
luminosity at fixed halo mass. At variance with respect to the be- 
havior of the UV magnitude, the Lya luminosity of hosts in larger 
halos is already declining for ages less than tait at 2 > 6, since Lya 



Long GRBs and Their Hosts at High z 7 




Figure 8. The comoving number density of galaxies as a function of UV magnitude at 1350 A Afisso (left) and of SFR (right). In botli panels solid, dashed 
and dotted lines refer to redshifts z = 3, 6, and 10, respectively. The sets of blue and green curves refer to no metallicity threshold and Zcrit = Zq/W. 
The data points reported in the left hand panel are from Steidel et al. (2001; circles) and Yoshida et al. (2006; squares) at 2 ~ 3 — 4, Yoshida et al. (2006; 
diamonds) at z ~ 4 — 5, and Bouwens et al. (2006; stars) aX z ~ 5 — 7. 



emission is sensitive not only to dust but also to neutral hydrogen 
absorption. 

Anyhow we expect a strict correlation between GRB hosts 
and Lyman Break Galaxies (LEG) and Lya emitters (LAE). As 
the galaxy model here adopted reproduces the high-redshift UV LF 
of LBGs and LAEs (see Mao et al. 2007), it is meaningful to esti- 
mate the UV LF expected for GRB hosts, namely by imposing the 
condition Z< Zcrit ~ 0.1 Zq. The result is shown in Fig. 8. In 
general the conditional UV LF depends on the ratio between tait 
and the time during which the host is brighter than a fixed abso- 
lute magnitude. As expected, the metallicity cut depresses the UV 
luminosity distribution, but this depression is minimal at the faint 
and bright ends, and maximum at intermediate values, as expected 
from the effects of the cut on i\fi35o in individual systems. 

In particular, the time spent by galaxies in large halos as LBGs 
- before turning into strongly submm emitting galaxies because of 
dust attenuation - is short and comparable with fcrit- At the faint 
end the luminosity of small hosts can fall below a fixed magnitude 
in a relatively short time, as shown in Fig. 8 for hosts in halos with 
Mh ~ 10^" Mq. It is also apparent that intermediate mass halos 
and the corresponding intermediate UV luminosity objects are sig- 
nificantly more affected by the condition on metallicity. Obviously 
at high redshift, e.g. z ~ 10, the cut reduces the visibility timescale 
only by a small fraction. 

We can conclude that GRB hosts should well reproduce the 
LBG LF at its faint and bright ends. Similar conclusions hold for 
LAE LR 

The shortness of tait ~ 50 Myr implies that the hosts must 
exhibit quite large [a/Fe] enhancement. In fact, the cumulative 
fraction of SNIa explosions after an instantaneous burst of star for- 
mation (1 at 12 Gyr) is negligibly small for 50 Myr (corresponding 
to the lifetime of a 7 A/0 star), for a wide class of progenitor mod- 
els, i.e. Single Degenerates (SD), and Double Degenerates (DD) 
exploders (Greggio 2005). At 0.1 Gyr it reaches 10 per cent in the 
most favorable case of close DD scheme. Furthermore, the typi- 
cal timescales required to significantly decrease the [a/Fe] ratio 
from its initial value produced by a single generation of massive 
stars (approximately 0.4 — 0.6 depending on metallicity, IMF and 
stellar yields, Gibson et al. 2003), is 0.3, 1 and 3 Gyr, for close 
DD, SD and WIDE DD progenitors, respectively. Thus while in the 



most massive ellipticals the duration of the burst of star formation 
(around 0.3 Gyr) may be enough for a mild pollution by Fe-peak 
elements (depending on the assumed scenario for SNIa progeni- 
tors), the lower timescale required to reach the critical metallicity 
implies that GRB hosts should display the original pattern of heavy 
elements produced by massive star chemical evolution. Notice that 
the subsequent history of star formation decreases the initial a- 
enhancement more in smaller than in larger objects (see Fig. 1). 
We conclude that at high redshift, GRB hosts (and the similarly 
young LAEs) exhibit the highest values of a-enhancement, even 
higher than those of high-z quasar hosts and of their descendents, 
the massive elliptical galaxies. 

Finally, we can also estimate the distribution of intrinsic SFRs. 
It is worth stressing that generally the hosts with SFRs exceeding 
a few X 10^ Mq yr~^ are significantly less than the overall galaxy 
population. This suppression reflects not only the steep halo mass 
function, but also the short duty cycle of GRB hosts following the 
rapid metal enrichment in the initial stages of star formation. And 
again the effect gets more pronounced at decreasing z due to the 
increased fraction of observable hosts which are below threshold. 
On the other hand, the model allows for the existence of GRB 
hosts with SFRs exceeding a few xlO^M© yr ~^ within halos 
endowed with Mh ^ lO^'^ Mq : as shown in Fig. 1 and Fig. 8, the 
fraction of these hosts should amount to a few per cent of the to- 
tal. Nevertheless they should not have formed a large amount of 
stars M* < 10^° Mq, exhibit specific SFR M*/M* > 10"* yr^^ 
and should be rather extincted by dust (see Fig. 1). 



4 DISCUSSION 

The adopted galaxy formation scenario coupled with the metal poor 
collapsar model suggested by stellar evolution have been exploited 
to infer the above results, which include GRB counts and redshift 
distribution and the complete description of the relevant properties 
of their hosts, such as SFR, mass in stars, chemical evolution in the 
cold star forming gas and in the stellar component within individual 
galaxy halos, with specified mass and formation redshift. In this 
§ 4 these results are discussed in the light of presently available 
observations and they are compared to previous studies. 
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4.1 GRB progenitor rates and GRB counts 

The imposed metallicity threshold affects the GRB progenitor red- 
shift distribution. Instead of peaking at around z « 3 as in the case 
without metallicity constrain, the threshold Zcrit = Zq/1Q yields 
a broad peak at 2 « 6 (Fig. 2). A significant fraction of progenitors 
are expected at high redshift: approximately 60 per cent at z > 4 
and approximately 30 per cent at « > 6 (see Fig. 3). This behav- 
ior is shared by all models which associate GRBs with metal poor 
progenitors (e.g. Natarajan et al. 2005; Langer & Norman 2006; 
Salvaterra & Chincarini 2007). 

In particular, our findings are similar to those inferred by 
Langer & Norman (2006), who explored the effect of a metallicity 
threshold in the environment of GRBs, by adopting a mass-stellar 
metallicity relation for galaxies and an average cosmic metallic- 
ity [Z] dependence on redshift d[.Z]/d2: ~ —0.15 dex. This law, 
derived to reproduce the metallicity of stars in galaxies, has been 
extrapolated to the ISM. This extrapolation implies that at 2: > 7 all 
cosmic star formation occurs in environments with [Z] < 0.1 Zq, 
independently of the processes occurring in galaxy halos. As a re- 
sult the corresponding rate of low metallicity SNae peaks at « « 6. 
In our model the combination of cooling and feedback processes 
implies that at 2 w 7 the SFR in cold gas with Z <Q.l Zq amounts 
to about 40 per cent of the cosmic one, due to the short timescale for 
gas enrichment. Therefore, despite of the quite different approach, 
the redshift distributions of GRB progenitors found by Langer & 
Norman (2006) and Yoon et al. (2006) are very similar to ours, 
since their assumed SFR for metal poor stars peaks at z « 6, simi- 
larly to our SFR in cold gas environment (see Fig. 2). 

The combination of the galaxy formation scenario with the 
low metallicity collapsar hypothesis for GRB events leads to a good 
agreement of the predicted and observed (bright, S>1 ph s~^ 
cm~^) 5W7Fr counts (see Fig. 4). This naturally follows - without 
any tuning - from the derived number of GRB progenitors, under 
quite simple and plausible assumptions on the GRB jet opening an- 
gle, prompt gamma-ray luminosity distribution and spectral shape, 
without requiring any luminosity evolution (see § 2.3 and § 3.2). 

The corresponding redshift distribution (see Fig. 5) implies 
that at this bright limit the adopted model predicts about 6 GRBs 
per year, while only one GRB at « > 6 has been identified in two 
years of SWIFT operation. As a matter of fact, only for about half 
of the SWIFT bursts an optical afterglow has been observed, and 
for only 30 per cent it has been possible to infer a redshift estimate. 
It is therefore quite reasonable to guess that this fraction could be 
about 10 per cent for bursts at z > 6 (see Fiore et al. 2007). 

Recently Salvaterra & Chincarini (2007) obtained a good for- 
mal fit to the SWIFT counts, down to S* ~ 0.4 ph s^^ cm~^, adopt- 
ing as free parameters the count normalization and the GRB LF 
(two further parameters). They considered the case of low metal- 
licity environment by adopting a kinematical model and predict the 
occurrence of only 1 GRB yr^^ for the SWIFT bright flux limit 
(1 ph s~^ cm~'^). Their result imply a redshift determination effi- 
ciency for GRBs greater than 50 per cent. 

It is worth noticing that reaching completeness down to 0.1 
ph~^ s~^ cm~^ would significantly increase the number of de- 
tected GRBs and in turn allow to explore the Universe during the 
recombination epoch, « > 8, with good statistical significance. A 
further decrease in the flux limit to 0.01 ph~^ s~^ cm~^, would 
instead only increase the GRB sample by a factor of 2; at this flux 
limit practically all GRBs would be detected (see Figs. 5 and 6). 

In order to account for the trend of GRBs to be at substan- 
tial redshift Firmani et al. (2005) proposed that the GRB LF is 



evolving. Daigne et al. (2006) tested the hypothesis of an increas- 
ing efficiency of GRB production by massive stars with increasing 
redshift. Detailed redshift distributions of GRBs will discriminate 
among these different possibilities. 

We can conclude that the hypothesis that metal-poor, rapidly 
rotating massive stars are the GRB progenitors (Woosley & Heger 
2006; Yoon et al. 2006) is consistent with the observed SWIFT 
counts. Clearly determinations of GRB redshifts will be extremely 
informative on the progenitor and galaxy formation models. 

We stress that the adopted galaxy formation scenario exploits 
quite a standard IMF, independent of the gas metallicity. As a mat- 
ter of fact, we showed that the cold star forming gas is rapidly {t <C 
10^ yr) enriched to the possible threshold around 3 — 5 x 10~* Zq , 
below which the IMF might be strongly biased toward high mass 
stars (see Bromm et al. 2001; Schneider et al. 2006). Therefore the 
possible contribution of pop III stars with IMF strongly biased to- 
ward high masses is not considered here. However, Bromm & Loeb 
(2006) showed that at z ~ 10 the contribution of pop III to cosmic 
SFR could be of order of 1/10 of the overall SFR, becoming domi- 
nant at z > 15; the GRB rate from pop III massive stars would grow 
correspondingly. 

4.2 Properties of Long GRB Host Galaxies 

Though the idea that GRBs are preferentially located in metal poor 

environments is attractive, nevertheless observational estimates of 
the metal content of host galaxies are still controversial. While most 
of the results suggest that at high-z GRB hosts exhibit metallicity 
Z < 0.1 - 0.3 Zq (Prochaska et al. 2007, Price et al. 2007), there 
are claims of higher metal content (see Savaglio et al. 2003). 

Once a metal poor (Z ^ 0.1 Zq) collapsar model is assumed, 
our galaxy formation scenario predicts that GRB hosts are very 
young, with age less than 5 x 10^ yr, independently of the halo 
mass. This young age directly mirrors the predicted short timescale, 
tcrit ~ 5 X lO'^ yr, of chemical enrichment of the cold gas, indepen- 
dently of mass. Such independence makes the GRB rate at high-z a 
good tracer of the virialization rate of relatively small galaxy halos, 
Mh < X 10" Mq (see Fig. 6). 

A definite prediction issuing from their youth is that GRB 
hosts should exhibit high [a/Fe] -enhancement, as their metal con- 
tent directly reflects the chemical yields of core collapse SNae. In- 
deed ,in the most favorable scenario of close DD progenitors, type 
la SNae may halve the [a/Fe] ratio produced by the generations of 
massive stars in a few 10* yr, somewhat longer than tcvit- Recent 
observations by Prochaska et al. (2007) suggest that a/Fe ratios are 
more than 3 times the solar value. Though differential depletion 
could be responsible for the result, the young age of the hosts is a 
much more palatable explanation. 

The shortness of tait implies that star formation has not 
much proceeded, i.e. the model predicts high specific star forma- 
tion Mi, /Mi, >2 X 10^*^ yr^\ Notice that even the specific star 
formation is almost independent of halo mass. Specific SFRs at 
this high level have been claimed by several authors (Fruchter et 
al. 1999; Le Floc'h et al. 2003; Christensen et al. 2004; Savaglio et 
al. 2006; Castro Ceron et al. 2006; Michalowski et al. 2008). We 
caution that metal poor gas may be left over particularly in external 
regions of relatively old galaxies, hosting a burst of star formation 
even if most stars already formed and the gas was already metal en- 
riched. However, we expect that relatively evolved hosts endowed 
with significant stellar mass are exceptions: GRB 020127 could be 
one of such cases (Berger et al. 2007). The vast majority of GRBs 
are hosted by small galaxies soon after their first stars shine. 
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Since the hosts are very young galaxies, our model predicts 
that they are gas rich objects with large column density A'h, once 
more independently of their galaxy halo mass. Observations con- 
firms the tendency for GRB hosts to exhibit large /Vh (Prochaska et 
al. 2007; Schady et al. 2007). At odd with these results, Tumlinson 
et al. (2007) set stringent upper limits on the molecular hydrogen 
(Hi) abundance, concluding that the fraction of molecular hydro- 
gen H2 /HI is extremely low in the examined hosts. These authors 
also point out that the deficiency may be related to low dust abun- 
dance, H2 formation being catalyzed on the surface of dust grains. 
An additional possibility is the destruction of H2 by UV radiation. 
The model here proposed predicts that GRB hosts are dust poor 
and are pervaded by a intense UV radiation field; therefore they are 
expected to be quite poor in molecular hydrogen. 

Concerning the UV emission, most of the hosts 
(Mh<1O^"M0) should have UV magnitude at 1350A in 
the range —20 < M1350 < — 16. The most luminous hosts are 
bright enough to be included in the already available LBG LF. 
Interestingly, Jakobsson et al. (2005) have shown that GRB hosts 
are tracing the faint end of the LBG LF at 2 « 3. 

The predicted host Lya luminosity should fall in the interval 
2 X 10*° < LLya ^ 2 X 10*^ erg s~^, only marginally overlapping 
with the range explored by currently available Lya LF (Fig. 7). Tan- 
vir & Levan (2008) have shown that the UV rest frame luminosity 
distribution of Lya selected galaxies and GRB hosts at 2; 3 are 
quite similar. We notice that Mao et al. (2007) pointed out that LAE 
are expected to be younger, with lower stellar masses, more com- 
pact and associated with less massive halos than LBGs. 

In summary, observations supports the model prediction that 
the GRB hosts trace the faint end of the LF of LBGs and LAEs. 

A further issue concerns the amount of dust in GRB hosts. 
While most hosts have not been detected at mid-IR or submm 
wavelengths (e.g. Tanvir et al. 2004; Le Floc'h et al 2006; Priddey 
et al. 2006) in a few cases submm and radio emission have been 
actually detected (Michalowski et al. 2008). How can the existence 
of these hosts be interpreted within the proposed scenario? 

The model predicts that only GRBs hosted in large galac- 
tic halos Mh > lO'^^ Mq can have significant dust absorption. 
These also exhibit large SFR (M* ^100 Mq yr~^) and relatively 
small stellar mass (M* < 10^° Mq). These properties correspond 
to those inferred by Michalowski et al. (2008) for the submm and 
radio detected objects. The model also predicts that these systems 
should represent only 1/20 of all the GRB hosts at 2 « 1 — 2 and 
practically disappear at a > 5. 

We stress that the conclusions regarding the host properties 
refer to the average GRB population, and do not exclude the pos- 
sibility that individual GRBs might reside in low metallicity local 
regions within their host (as might be the case for GRB 060206, 
Fynbo et al. 2006). 

The above findings have been derived in the framework of a 
galaxy formation scenario, where a key role is played by the energy 
feedbacks provided by SNae and quasars. However, the fast chem- 
ical enrichment leading to gas metallicity above Zcr\t ~ 0.1 Zq 
takes place at early galactic times, before quasar feedback becomes 
effective. On the other hand, the stellar feedback is very relevant 
since it regulates the star formation activity even at early galactic 
times. 



5 SUMMARY 

We have explored the cosmological consequences of the assump- 
tion that metal poor and rapidly rotating single stars are the progen- 
itors of most long GRBs. Our main conclusions are: 

• The overall long GRB rate amounts to approximately 300 
yr~^ sr~^. Bright SWIFT counts are reproduced by assuming a 
non-evolving prompt (gamma-ray) luminosity distribution. 

• Above a flux limit of 1 ph s^^ cm~^ about 30 per cent of 
GRBs are predicted to be at 2 > 6 and 10 per cent at 2 > 8, amount- 
ing to approximately 13 for two years of SWIFT operation. Only 
one have been located above 2 > 6 in two years: this would require 
a redshift determination efficiency around 10 per cent, to be com- 
pared with the 30 per cent efficiency for events at lower z. 

• The host galaxies are very young, with age less than 5 x lO'^ 
yr, gas rich (large column densities) but poorly extincted systems 
(^v ^ 0.1), because of their chemical immaturity. Only the few per 
cent of hosts associated with large halos (Mh lO^'^ Mr.;) have 

large extinction (^v > 0.3), high SFR (M* > 100 Mq yr"^) and 
can be detected at relatively bright submm flux levels. This result 
has impUcations for the origin of 'dark' GRBs [about 1/3 of SWIFT 
GRBs can be considered 'dark', e.g. Schady et al. (2007)], lack- 
ing a detection of the optical afterglow. Dark GRBs should largely 
comprise a population of high 2 events, rather than highly extincted 
systems. 

• The young age of hosts implies that (i) the specific SFR is 
high (M*/M* > 2 X 10"*) and (ii) the ratios of abundance of dif- 
ferent chemical elements are just those of the respective chemical 
yields of type II SNae, i.e. large a-enhancements should be the rule. 

• Most of the hosts (10^ < Afn < lO" Mq) have UV mag- 
nitude in the range —20 < M1350 < — 16 and Lya luminosity in 
the range 2 x 10"° < Li^ya < 2 x lO"^ erg s"\ They trace the for- 
mation of small galaxies in small halos, and as a consequence the 
faint end of the LBG and LAE LF. These hosts would reionize the 
Universe at 2 « 7 (see Mao et al. 2007). 
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APPENDIX A: A SIMPLE RECIPE FOR STAR 
FORMATION AND DUST OBSCURATION IN 
PROTOGALAXIES 

Granato et al. (2004) have proposed a model for early galaxy for- 
mation, in which the most relevant processes (gas cooling and in- 
flow, star formation and gas accretion onto BH, stellar and quasar 
feedback, gas outflow) are included in a set of equations, that can 
be solved with straightforward numerical computations. In this Ap- 
pendix we present a simplified version of the model describing the 
SFR, mass in stars and chemical evolution, which are relevant for 
this work. The analytical formulae presented below are very good 
approximations of the results found by solving the full set of equa- 
tions of the model (for details see Lapi et al. 2006 and Mao et al. 
2007). We stress that, because of their fundamental character, the 
equations listed below describes the main aspects of star formation 
and chemical evolution in protogalaxies at high redshift. 

When a DM halo of mass Mh reaches the virial equilibritun, 
it contains a mass Minf (0) = /cosm A/h of hot gas at the virial 
temperature, /cosm « 0.18 being the mean cosmic baryon to DM 
mass-density ratio. The gas in virial equilibrium flows toward the 
central region at a rate Mcond = Mnf/^cond where the conden- 
sation timescale tcond = max[tcooi(i?H), tdyn(i?H)], is the maxi- 
mum between the dynamical time and the cooling time at the halo 
virial radius i?H. When computing the cooling time, a clumping 
factor C in the baryonic component is also included: C7 > 7 im- 
plies tcooi(i?H) "^idyniRii) OH relevant galaxy scales at high- 2:. 
By defining such a condensation time, we implicitly neglect the ef- 
fect of angular momentum. However, angular momentum decays 
on a dynamical friction timescale tu¥ « 0.2 (5/ in 5) idyn, where 
^ = Mh /Mc and Mc is the typical mass cloud involved in major 
mergers (e.g. Mo & Mao 2004); major mergers, which are very fre- 
quent at high redshift and in the central regions of halos, imply ^ ~ 
a few. 

The model also assumes that quasar activity removes the hot 
gas from the halo through winds at a rate °; the equation for 
the diffuse hot gas is then 

Minf = -Mcond - MiSf ° . (Al) 

The cold gas is piled up following the cooling of hot gas, is 
consumed by star formation (M*), and is removed by the energy 
feedback from SNae (,Mco\d) and quasar activity (M^^fd°): 

Mcoid = Mcond - (1 - n)M^ - McoS - M^oiA - (A2) 
where TZ is the fraction of gas restituted to the cold component by 
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the evolved stars. Under the assumption of instantaneous recycling, 
7^ f» 0.54 for the adopted IMF (7^ 0.3 for a Salpeter IMF); this 
value of TZ is an upper Umit, since only a fraction of evolved stars 
have a significant mass loss in the evolutionary phases considered 
here. However, the relevant results are only very weakly sensitive 
to the chosen value, in the physically allowed range. The mass of 
cold baryons that is going to be accreted onto the central supermas- 
sive BH is small enough to be neglected in the above equation (see 
Granato et al. 2004). 

Stars are formed at a rate 

dMcoid Mcoid 



(A3) 



max[tcool,tdyn] U 

where now tcooi and tdyn refer to a mass shell dMcoia, and U is 
the star formation timescale averaged over the mass distribution. 
The rate of cold gas removal due to SNae is parameterized as 



(A4) 



where the averaged efficiency 

^^SN €SN £^SN 



/3s N 



V0.05/ 



0.6 



Mu 



1.4 X 10-VMq 



-2/3 



7 



(A5) 



depends on the number of SNae per unit solar mass of condensed 
stars nsN, the energy per SN available to remove the cold gas 
esN £^SN, and the specific binding energy of the gas within the 
DM halo, Ebind- Following Zhao et al. (2003) and Mo & Mao 
(2004), the latter quantity has been estimated, for z > 1, as i?bind = 

Vi f{c) (1 + /cosm)/2 ~ 5.6 X 10^* (Mh/W" Msf/'-' [(1 + 
z)/7] cm^ s^^. Here Vh is the halo circular velocity at the virial 
radius and /(c) ^ 2/3 + (0/21.5)"'^ ~ 1 is a weak function of the 
halo concentration c ^ a few. Lapi et al. (2006) have shown that 
high redshift LFs of quasars and galaxies constrain esN w 0.05; 
the same value is required in order to reproduce the fundamental 
correlations between local ellipticals and dormant BHs. 

By analyzing the results of the numerical solution of the full 
set of equations by Granato et al. (2004), it is apparent that the 
term of quasar feedback is important only during the final stage of 
BH growth, around 2 — 3 e— folding times (approximately 10* yr) 
before the peak of quasar luminosity, when the energy discharged 
by the quasar is so powerful to unbind most of the residual gas, 
quenching both star formation and further accretion onto the super- 
massive BH. The time integral over the quasar bolometric power 
exceeds the gas binding energy after 

Atbnrst « 2.5 X 10* ( T ( ,,5". ^ yr > (A6) 



7 



where J-{x) = 1 for a; ^ 1 and J-{x) 



10^2 Afg 

for a;- < 1. Therefore 
a good approximation for the star formation history is obtained by 
neglecting the quasar feedback effect in Eqs. (Al) and (A2) and by 
abruptly stopping star formation and accretion onto the central BH 
afterAtburst since halo virialization. 

Then Eqs. (Al) and (A2) can be easily solved, with the out- 
come that the infalling mass declines exponentially as Minf (t) = 
Mnf (0) e~*/*'=°"'^, while the SFR evolves according to 

Minf(O) 



lis) i 



Q */*cond g 



-d] , (A7) 



tcond(7 - 

with 7 s 1 — 7i -I- /3sN. The quantity s = tcond/** is the ratio be- 
tween the timescale for the large-scale infaU estimated at the virial 



radius and the star formation timescale in the central region; it cor- 
responds to s ~ 5, both for an isothermal or NFW (Navarro et al. 
1997) density profile. 

The following expression well approximates the condensation 
timescale (see Mao et al. 2007): 



tcond « 4 X 10' 



l + z 
7 



1012 Mq 



yr . (A8) 



The scaling with redshift reflects the dependence of the dynami- 
cal time; the weak dependence on Mh reproduces the impact of 
the energy feedback from the quasar on the infalling gas, which is 
stronger for more massive halos hosting more massive BH. 

In order to compute the metal content of the cold gas, one has 
to take into account the infall of primordial abundance gas, the en- 
richment due to earlier generations of stars, and outflows driven by 
winds generated by SNae and quasars. We assume, as common, that 
in the cold gas there is an instantaneous mixing of metals released 
by stars. The corresponding simplified equation reads 

Zeo.d(t) = ^4^#^ + ¥^, (A9) 



^cond M;old {t) t* 



where 



dm* m* qz {rrii,) <j){mi,) 



(AlO) 

The IMF is denoted by 4>{mi,) and qz{mn) is the metal yield of 
stars of mass m*. From Fig. 1 it is apparent that metallicity rapidly 
(in a time less than a few 10* Gyr) increases from primordial con- 
tent to the limiting value 



S7 ■ 



-Az, 



(All) 



with Az = 0.043 for the adopted IMF (Az = 0.021 for the 
Salpeter IMF). 

In early galaxy evolution an important role is played by dust, 
that absorbs the UV emission and re-radiates it in the mid and far- 
IR band. The amount of dust in a galaxy is expected to be correlated 
with that of cold gas [or SFR, see Eq. (A3)] and with metallicity. 
Mao et al. (2007) have shown that it is possible to describe the 
luminosity-reddening relation found by Shapley et al. (2001) for 
a « 3 LBGs with the simple law 



0.35 



Moyr- 



Z 



(A12) 



where .4i:ir,() is the attenuation at 1350 A. By inserting in the above 
equation the SFR and metallicity, we get the attenuation as function 
of time. Then the UV magnitude is simply given by 



Afi350 ~ -18.6 - 2.5 log I 



Mq yr-\ 



+ Al350 . (A13) 



This simplified treatment of dust attenuation proved to be 
quite a good approximation for low-luminosity LBGs and LAEs, 
which exhibit low attenuation (Mao et al. 2007). 

On the other hand, for massive halos attenuation is large soon 
after 10* yr; this corresponds to the quick appearance of very lu- 
minous submm-selected galaxies. For these systems the model by 
Granato et al. (2004) includes a sophisticated treatment of dust at- 
tenuation through the GRASIL code (Silva et al. 1998). 

The scheme presented here traces the evolution of single 
galaxy halos as function of time, given their mass and their for- 
mation redshift. The results can be interfaced to the formation rate. 
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d? NsT /dM-a At^ (see § 2.2) yielding LFs, counts and redshift dis- 
tribution for spheroidal galaxies. Moreover, the quasar LF can also 
be reproduced (Lapi et al. 2006), as well as the LF of LBGs and 
Lya emitters (Mao et al. 2007). A summary of the parameters used 
by the model and of its achievements is presented in Tables 1 and 
2ofLapietal. (2006). 



